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ABSTRACT 

We explore, through a simplified, semi-analytic model, the formation of dense 
clusters containing massive stars. The parent cloud spawning the cluster is rep- 
resented as an isothermal sphere. This sphere is in near force balance between 
self-gravity and turbulent pressure. Self-gravity, mediated by turbulent dissipa- 
tion, drives slow contraction of the cloud, eventually leading to a sharp central 
spike in density and the onset of dynamical instability. We suggest that, in a real 
cloud, this transition marks the late and rapid production of massive stars. 

We also offer an empirical prescription, akin to the Schmidt law, for low- 
mass star formation in our contracting cloud. Applying this prescription to the 
Orion Nebula Cluster, we are able to reproduce the accelerating star formation 
previously inferred from the distribution of member stars in the HR diagram. 
The cloud turns about 10 percent of its mass into low-mass stars before becoming 
dynamically unstable. Over a cloud free-fall time, this figure drops to 1 percent, 
consistent with the overall star formation efficiency of molecular clouds in the 
Galaxy. 

Subject headings: open clusters and associations: individual (Orion Nebula Clus- 
ter) — stars: formation — stars: pre-main-sequence — ISM: clouds 



1. Introduction 

There is growing evidence that the formation of stellar groups is a relatively slow process. 
More specifically, a star cluster appears within its parent molecular cloud over a period long 



compared to the cloud's free-fall time, as gauged by the mean gas density. iTan et al.l (120061 ) 
have summarized several lines of argument leading to this conclusion. The gas clumps 
believed to form massive clusters look round, indicating that they are in force balance, and 
not a state of collapse. Massive clusters themselves have smooth density profiles, again 
in contrast to a dynamical formation scenario. The observed flux in protostellar outflows 
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indicates a slow accretion rate, and therefore a long star formation time scale. Finally, 
placement of young clusters in the HR diagram yields age spreads in excess of typical free- 
fall times (see also Palla & Stahler 2000). 

Many researchers have performed direct numerical simulations of molecular clouds; their 
results also bear on the issue of the star formation time scale. In a typical simulation, the 
computational volume is filled with a magnetized, self-gravitating gas that has a turbulent 
velocity field. If the turbulence is only impressed initially, it dies awa y in a crossing time , 
and most of the gas condenses into unresolvably small structures (e.g., iKlessen et al.lll998l ). 
Since the crossing and free-fall times are similar in a molecular cloud, so me authors have 



concluded that all clouds produce stars rapidly, while in a state of collapse (|Hartmann et al. 



2OOII ). Others have used empirical arguments to make the same point (lElmegreeru I2OOOI ) . 
This view is at odds with the observations concerning cluster-forming clouds cited above. 
Moreover, the simulations show that, if turbulence is driven thr oughout the calculation, t he 
rate of star formation can be reduced to a more modest level ( Mac Low fc Klessen 2004j )F] 
It is plausible that the turbulence is indeed driven by the cloud's self-gravity, a point we 
shall amplify later. 

The emerging picture, then, is that molecular clouds both evolve and create internal 
clusters in a quasi-static fashion. That is, the structure as a whole is nearly in force balance, 
until it is eventually destroyed by the ionizing radiation and winds from the very stars 
it spawns. The inferred masses of all clouds larger than dense cores greatly exceed the 
corresponding Jeans value, evaluated using the gas kinetic temperature. Thus, self-gravity 
must be opposed by some force beyond the relatively weak thermal pressure gradient. The 
extra support is generally attributed to MHD waves generated by internal, turbulent motion 
( lElmegreen fc Scald |2004| ). This motion, which is modeled in the numerical simulations 
just described, impr ints itself on molecular line transitions, giving them their observed, 
superthermal width (lArons fc MaxI 1 19751 : iFalgarone et al.l Il992l ) . 



In this paper, we follow the quasi-static contraction of a spherical, cluster-forming cloud. 
Contraction is facilitated by the turbulent dissipation of energy. This investigation continues 
and extends an earlie r one that was part of our study of the Orion Nebula Cluster (ONC) 
(iHuff fc Stahlerl |2006| . Paper I). Here we track in more detail the changing structure of a 
generic cloud, taken to be in near balance between self-gravity and turbulent pressure. We 
find that contraction eventually causes the density profile to develop a sharp, central spike. 



^The actual rate of condensation depends on the mag nitude of the sonic length, i.e., the size scale of turbu- 
lent eddies whose velocity matches the local sound speed (| Vazquez- Semadeni et al 1200,4 iKrumholz fc McKe^ 
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Suc h a region of growin g; density is a plausible environment for the birth of massive stars 



see 



Stabler et al.ll2000f ) 



We also track, using a simple, empirical prescription, the formation of low-mass stars 
in our contracting cloud. For reasonable parameter values, stellar births occur throughout 
the cloud over a period of order 10^ yr. The global rate of star formation rises with time 
monotonically, i.e., the formation accelerates. Extended, accelerating pro duction of stars 



is als o found empirically when one analyzes clusters in the HR diagram (jPalla fc Stabler 



20001 ) . Indeed, it is not difficult to match specifically the global acceleration documented in 
the ONC. Here, the star formation rate depends on cloud density in the same manner as the 
classic Schmidt law. 

In Section 2 below, we present the basic physical assumptions underlying our model. 
We also give a convenient, nondimensional scheme. In Section 3, we introduce our treatment 
of turbulent dissipation, and calculate the interior evolution of the cloud as it contracts 
toward the high-density state. Section 4 offers our prescription for low-mass star formation, 
and compares the resulting birthrate with the ONC data. Finally, Section 5 discusses the 
broader implications of our findings, as well as their utility for future work. 



2. Formulation of the Problem 



2.1. Physical Assumptions 



We focus on molecular cloud clumps that are destined t o produce the highes t-density 
clusters, i.e., those containing massive stars near their centers. IShirley et al.l (120031 ) used CS 
observations to study a sample of 63 clumps already containing massive stars, as evidenced 
by water maser emission. These clouds are nearly round, with median projected axis ratios 
of 1.2. It is thus a reasonable approximation, and certainly a computationally advantageous 
one, to take our model cloud to be spherically symmetric. 



The clumps observed by lShirley et al.l (120031 ) have a median radius of 0.32 pc and mass 
of 920 Mq. a cloud of this size and mass has internal, turbulent motion well in excess of 
the sound spe ed, where the latter is based on the typical gas kinetic temperature of 10 K 
(jLarsonlll98ll ). This bulk motion excites a spectrur n of MHD waves, i.e., pe rturbations to 
the interstellar magnetic field threading the cloud (IFalgarone fc Pugetl 119861 ). Such waves 
exert an effective press ure that can, at least in principle, support the cloud against global 
collapse JPudritzlll99"o[ ). 



Fatuzzo &: Adams 



(119931 ) studied the mechanical forcing due to MHD waves propagat- 
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ing in a one-dimensional, self-gravitating slab. They considered two slab with an 

embedded magnetic field oriented parallel to the slab plane, and one with an internal field in 
the normal direction. In the first case, Fatuzzo & Adams showed that magnetosonic waves 
provide a normal force. In the second, it is Alfven waves that exert the force, also in the 
normal direction. Thus, Fattuzzo & Adams verified explicitly that the waves counteract 
gravity, even in the absence of wave damping. 



McKee fc Zweibell (119951 ) extended this result. Using the pioneering analysis of iDewar 
(Il970l ). they showed that Alfven waves generated by a turbulent wave field exert an isotropic 
pressure, regardless of the background geometry. McKee & Zweibel derived a simple depen- 
dence of the wave pressure P on the local density: 

P oc _ (1) 

If we ignore the relatively small thermal pressure, then the cloud can be described as an 
n = —2 polytrope. 

Are the structures of real clouds consistent with this polytropic wave pressure? One, 
indirect, argument indicates they are not. McKee & Zweibel also demonstrated that P is 
proportional to p times the square of the (randomly oriented) velocity fiuctuation 6v. It 
follows that 

5v oc p-^/^ (2) 

in this model. Now the density in an n = —2 polytrope tends to approach a power law 
outside the central plateau, such that p is proportional to r~^/^. From equation (2), it 
follows that 6v is proportional to r^/^. 

Consider the nearly spherical cloud, now gone, that produced the ONC. This cloud 
was recently driven of f by the Trapezium stars, which themselves have ages of about 10^ yr 



( iPalla fc Stabler! I2OOII ). The disruption itself occurred well within the cluster crossing time 
of about 10^ yr. Hence, the present-day velocity dispersion of the stars should refiect the 
prior Sv of the gas. But the dispersion of ONC proper motions has negligible variation from 



the center to the outskirts of the cluster (j Jones fc Walked 1 19881 ). These measurements span 



at least a decade in radius, over which Sv should vary by a factor of 2.2, according to the 
polytropic relation. 

Our conclusion, based on this admittedly limited evidence, is that a more realistic 
model of the internal turbulence has a spatially constant velocity dispersionj^ If we further 



^Inside giant cloud comple xes, the observed velocity dispersion increases with the size of the substructure 
(jOssenkopf fc Mac Lowll2002f ). Again, we are focusing on a single clump, where such considerations do not 
apply. 
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appropriate the relationship between P and 5v derived by McKee & Zweibel, we are then 
positing an isothermal equation of state: 

P = pal (3) 

Here ot, the effective isothermal sound speed, is taken to be a fixed constant at a 
given instant of time. This same quantity varies temporally; indeed, this latter variation 
essentially drives the cloud's evolution. We emphasize that ax does not, as in ordinary 
gas dynamics, give the magnitude of random, microsopic velocities. Instead, this quantity 
represents, however crudely, the bulk motion of turbulent eddies; these eddies create the 
pressure P via MHD waves. 

Since we are modeling the cloud as an isothermal sphere, we face the familiar difficulty 
that its mass is infinite unless the configuration is bounded externally. We therefore picture 
the cloud as being surrounded by a low-density, high-temperature medium with an associated 
pressure Po- This latter quantity is also the pressure at the boundary of our spherical cloud. 
The cloud density at the boundary, po, is found from equation (3), given knowledge of a|,. 



2.2. Nondimensional Scheme 

The mathematical description of a self-g ravitating, isothermal cloud in hydrostatic bal- 



ance is well known (see lStahler fc Pallall2004j . Chap. 9). All structural properties follow from 



the isothermal Lane-Emden equation: 



=expM'). (4) 



with boundary conditions ip{0) = ip'{0) = 0. Here, ip is the dimensionless form of the 
gravitational potential (pg-. 

ifj = (f)g/al . (5) 

The nondimensional radius ^ is obtained from the dimensional r using G, a^, and the central 
density pc- 

Equation (4) was derived using both Poisson's equation and the condition of hydrostatic 
equilibrium. The latter may be recast as a relation between the density at any radius, p, its 
central value, pc, and the potential: 

p = pc exp {-ip) . (7) 
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The full, dimensional mass Mo follows by integration of p over mass shells, Using equation 
(4) to evaluate the integral, one finds 

M, = ^^(et\ . (8) 



Here, the subscript denotes the cloud boundary. Similarly, we shall use Ro for the radius of 
that point, where the internal cloud pressure has fallen to Pq- 

In this standard formulation, all nondimensional variables are defined through the basic 
quantities a^, pc, and G. Although the standard variables will remain useful, the scheme 
itself is not well suited to describing cloud evolution at fixed mass. For this purpose, we 
shall also utilize a second nondimensional scheme, based on Mq, Pq, and G. 

Let A be the new nondimensional radius, and a the nondimensional effective sound 
speed. These are defined through 

A = ° ... r , (9) 



and 



Q,2 = (^Q) 



We further denote as 5 the nondimensional density: 

G^I^Ml^ 

^ = p3/4 P ■ (11) 

Since we will be discussing temporal evolution, we define a nondimensional time through 



Mo 



1/4 



t . (12) 



It will be useful to relate new nondimensional quantities to old ones. Thus, equation 
(6) tells us A as a function of ^: 

For the central density appearing here, 6c, we use equation (7), evaluated at the cloud 
boundary: 

exp (V>o) 

4 = . (14) 

Finally, a itself may be written in terms of standard variables by using equation (8): 

= V4^(^^2^^ exp(^o/2) . (15) 
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3. Cloud Evolution 

3.1. Internal Structure 

We now consider a sequence of isothermal spheres of fixed mass, all embedded in the 
same external pressure. We may describe each structure using the new, nondimensional 
variables. The sequence is characterized by a single parameter, the center-to-edge density 
contrast; we shall denote this ratio as (3. From equation (7), (3 can also be written as 

13 = exp (V'o) (16) 

Since V'lO is a known function, there is a one-to-one correspondence between our fundamen- 
tal parameter (3 and .^o, the old, nondimensional radius. The potential if) increases monoton- 
ically with ^, so (3 likewise increases with ^o- The lowest value of j3 is unity, corresponding 
to = 0. 

The internal velocity dispersion a varies along our sequence. We may track this change 
through equation (15). Thus, for each selected /?, we first find ipo from equation (16). From 
knowledge of the function ^'(O) the corresponding ^o, as well as {dip jd^)^. Equation 

(15) then yields a. 

It is equally straightforward to obtain the internal density profile, 5(A), of any model. 
Knowing ipo and a, equation (14) gives the central density, 5c. Proceeding outward, equa- 
tion (13) gives the value of ^ corresponding to each A. Again using ^'(O' equation (7) 
yields the density ratio, 5 /5c = exp^ip). When we get to the edge, A = Ao, we find that 
5o/5c = [3~^ = exp{—ipo), in agreement with equation (16). 

Figure 1 displays graphically the change of the cloud's structure as a function of p. 
Here we have plotted the radius, A(/?), of selected mass shells. As expected, a shell in the 
deep interior monotonically shrinks. Other shells, however, turn around. With rising f3, an 
increasing fraction of the cloud mass starts to expand. Such expansion costs energy. Thus, 
configurations of very high (3 are not physically accessible, as we shall see. 



3.2. Enthalpy and Dynamical Stability 



The lower dashed, horizontal line in Figure 1 corresponds to a /3- value of 14.1. This is 
the Bonnor-Ebert den sity contrast . In the stand ard analysis, clouds of higher contrast are 
dynamically unstable (jEbertlll955l : lBonnorlll956l ). We recall, however, that this instability 
arises from perturbations of a cloud held at fixed temperature. In contrast, our sequence 
has varying effective sound speed. The Bonner-Ebert contrast no longer marks a stability 
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transition. Nevertheless, this value, which we denote as /3min, is still of interest. It signifies, 
at least in an approximate way, the point where self-gravity starts to overwhelm external 
pressure as the main compressive force. Our description of cloud evolution will henceforth 
focus on such gravity-dominated configurations, i.e., those for which /3 > Pram- 

To analyze stability in the present sequence, we first need to invoke thermodynamics. 
We showed in Paper I that energy dissipation in an isothermal cloud results in a decrease of 
the total enthalpy. Returning to dimensional notation, equation (AlO) stated 

where L is the luminosity. The enthalpy H is the generalization, to a self-gravitating gas, of 
the classic definition: 

H = £'therm + -Egrav + PoV , (18) 

where £'therm and E'grav are the thermal and gravitational potential energies, respectively, 
and V is the cloud volume. 

To evaluate £'therm, we employ the general expression for a nonrelativistic gas, (3/2) fP dV . 
Using equation (3) for P, this integral becomes (3/2)Moa|,. Instead of evaluating -Egrav di- 
rectly, we invoke the virial theorem, in the form 

-E'grav — — 2£^therm + ^ PoV . (19) 

After expressing the cloud volume in terms of the radius, the enthalpy is 

H = J-Moal+^-^PoRl. (20) 
If we define a nondimensional enthalpy h through 

then equation (20) becomes 

3 9 16 TT , , 

Figure 2 shows h along our sequence of clouds. Again, wc restrict ourselves to gravity- 
dominated configurations, for which /? > /3min. We also recall that {3 increases monotonically 
along the sequence. Plotted here against a, the enthalpy dips to a minimum, then spirals 
inward toward a point. The latter corresponds to the singular isothermal sphere. For this 
special configuration, it may be shown that — (7r/2)^/^ and Ao = (l/87r)^/^. Thus, the 
hmiting value of h is —0.187. 
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However, this limiting value is never reached in the course of evolution. As long as the 
cloud releases energy into space, so that L > 0, equation (17) tells us that the enthalpy 
declines. Thus, the last accessible configuration coincides with the minimum-enthalpy point 
in Figure 2. Numerically, we find that h = —0.50 for this cloud. The corresponding density 
contrast (3 is 370. 

Consider now two configurations with identical values of h, very slightly above the 
minimum. These clouds, like all those in the sequence, have the same mass. We may view 
them as extremal states attained by the minimum-enthalpy cloud in the course of a normal 
mode of oscillation. Here, we are assuming that the cloud radiates negligible energy during 
an oscillation period, so that h remains constant. The two endstates are in precise force 
balance; intermediate states depart only slightly from this condition. In the small-amplitude 
limit, the oscillation has zero frequency, and the unperturbed, minimum-enthalpy, state 
represents a stability transition. 

In summary, an isothermal cloud becomes dynamically unstable at a density contrast f3 
of 370, provided the global enthalpy i s held fixed duri ng any oscillatory perturbation. This 



important fact was first discovered by IChavanid (120031 ) in the course of a general analysis of 



isothermal configurations. § Note again the marked contrast with the traditional, Bonnor- 
Ebert result. The much lower critical density contrast in that case [jS = 14.1) arises because 
the cloud releases — and draws in — as much energy as necessary to remain isothermal, even 
during a single oscillation period. This assumption would be inconsistent with our picture 
that the cloud is quasi-statically contracting due to relatively slow, turbulent dissipation. 

The minimum-enthalpy state thus marks the natural endpoint of the cloud's evolution. 
We denote as /3max the corresponding density contrast, and display this limit as the upper 
dashed, horizontal line in Figure 1. Clouds with higher density contrast, including the 
singular isothermal sphere, are inaccessible. 



3.3. Turbulent Dissipation 

Although we have drawn a number of conclusions regarding the changing structure of 
our model cloud, we have yet to discuss its temporal evolution. The quasi-static contraction 
envisioned here is facilitated by the release of energy. This emission must arise at the shock 
interface between colliding, turbulent eddies. Typical fiuid speeds are the virial value, i.e.. 



■^Chavanis finds a slightly higher critical (3 of 390. His minimum enthalpy value, in our units, is 
h = -0.493. 



- 10 - 



less than 10 km s~^ for the clouds of interest. Hence, the shocks radiate through far-infrared 
and submillimeter photons from low-lying transitions of molecules. The cloud is optically 
thin to such photons. The luminosity L in equation (17) is thus generated from the full 
interior. 

Consider, then, a representative volume of the cloud. The numerical simulations men- 
tioned previously have modeled the dynamics of a magnetized gas subject to an impressed, 
turbulent velocity field. Even if the fluid disturbances begin as incompressible Alfven waves, 
effici ent mode conve rsion produces compressible MHD waves that steepen and shock (see, 
e.g.. lGoldsteinlll978l ). Signficant en ergy is dissipate d during the characteristic crossing time 



of the largest eddies. For example, iMac Lowl (119991 ) found that 
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'^turb 

A 



(23) 



Here e is the energy loss rate per unit mass of gas, V^urb the average (rms) eddy speed, 
and A the dominate wavelength of the impressed turbulence. The empirical constant rj was 
measured by Mac Low to be about 0.4. 

In the simulations, turbulence is impressed on an arbitrary scale . Indeed, the question 
of what drives the turbulence remains controversial (IMac Lowll2004l ). Here, we recall the 
key fact that the mean, interio r velocities match the virial value over a large range of cloud 
sizes and masses ( lLarsonlll98ll ). It is likely, therefore, that self-gravity constitutes the basic 
driving mechanism, although a quantitative model is still lacking. If this basic idea is correct, 
then the "dominant wavelength" in equation (23) should be comparable to the cloud size. 
We therefore adopt, as our expression for the cloud luminosity L, a mass- integrated version 
of this relation: 



7] 



(24) 



where the prefactor rj does not necessarily have the value found by Mac Low. The factor of 
4 in the denominator reflect s the fact that the largest mode corresponds to overall expansion 
or contraction of the cloud (IMatznerl 120021 ) . 



3.4. Quasi-Static Contraction 

To follow the cloud evolution in time, we use our prescribed luminosity to alter the global 
enthalpy. Combining equations (17), (20), and (24), we recast the result into nondimensional 
terms: 

da^ , 9 dAo , x 
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The dependent variables a and Ao are already known implicitly in terms of ^o- (Recall 
equations (13) — (15).) We may thus regard equation (25) as giving the dependence of r on 
this same quantity. Since r/ is still unknown, we use instead the combination t^t: 

d{r]T) _ 6Ao da^ _ 64 7rA^ d\o .2g^ 
d^o «^ d^o d^o 

We integrate this equation numerically, setting rjr = at = 6.5, the value at the Bonnor- 
Ebert density contrast, and ending at = 25, the minimum-enthalpy state. Over this 
evolutionary span, rjr increases by 0.96. 

Figure 3 shows that neither the effective sound speed nor the cloud radius vary greatly 
during this period. The former increases monotonically, with a fractional change of 10 percent 
by the end. The radius gently decreases most of the time. Just before the unstable state is 
reached, the cloud surface begins to swell, in agreement with Figure 1. 

The temporal change of the central density is much more dramatic. As seen in Figure 4, 
Sc increases slowly at first, and then accelerates strongly at the end. Given the behavior of 
the cloud radius, this rapid compression evidently involves a small fraction of the total cloud 
volume. The left panel of Figure 5 displays the evolution of the full density profile. From 
bottom to top, the associated values of rjT are 0, 0.15, 0.23, and the final value, rjr = 0.24. 

Because the inner portion of the cloud undergoes such rapid compression, one may 
question the basic assumption of quasi-static behavior. Are all mass shells really moving at 
subsonic velocity? At the very center, the velocity must fall to zero at each time. But what 
of shells just outside the center? 

Figure 6 shows the velocity profiles for the final three times depicted in the left panel 
of Figure 5. Since the temporal variable in equation (26) is t]t, we plot the ratio v/rja. 
For rjT = 0.60 and 0.92, the velocity is subsonic throughout, even if rj were unity. In fact, 
plausible //-values are less than unity, as indicated earlier; we shall see in Section 4 below that 
rj ^ 0.3 best matches the ONC data. In the final profile, corresponding to the minimum- 
enthalpy state, the peak velocities occur in the deep interior, and are mildly supersonic 
{v > a). 

4. Modeling the ONC 
4.1. Prescription for Star Formation 



At present, we have scant knowledge of how the birthrate of stars scales with t he prop- 
erties of the cloud medium spawning these objects. Within the solar neighborhood, ISchmidt 
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( 119591 ) found the star formation rate to be proportional to the square of the local density; 
here, the rate is measure per unit volume. Schmidt's law thus states that the formation rate 
per gas mass scales linearly with the cloud densityB 

On the scale of galactic disks, it is established that the formation rate per unit disk 
area rises as S", where S is the total gas surface density, and the exponent n is about 1.4 



( iKennicuttI Il998l ) . However, it is not straighforward to relate this important finding to the 
present study. Each areal patch in the galactic observations comprises numerous molecular 
complexes, any one of which is far larger than the clumps of direct interest here. 

Returning to our model, if the clump indeed undergoes slow contraction, then its over- 
all star formation rate must increase with time. It is plausible that the local rate within 
each mass shell rises with that shell's density p. Following Schmidt, we posit a power-law 
dependence: 

^* = f f-V • (27) 
h \PiJ 

Here, Th^, is the mass in stars forming per unit time, per unit cloud mass. The fiducial density 
and time, pi and ti, are those from equations (11) and (12), respectively: 

p3/4 

Pi = — 7^ — m (28a) 

t, = . (28b) 

Finally, the exponent n is to be set by matching to observations. We stress that the 
prescription in equation (27) applies to low-mass stars. The formation of massive objects is 
a separate phenomenon. In our model, this occurs only in the high-density, central region of 
the final, minimum-enthalpy state. 

Equation (27) contains a nondimensional efficiency factor e. Since only a fraction of the 
cloud mass turns into stars, we expect e to be well under unity. If, for simplicity, we assume 
this parameter to be the same in all mass shells, then integration of equation (27) yields the 
total mass per unit time in new stars: 

M* = P^ rp^^'r'dr. (29) 

*lPl ^0 



"^In Schmidt's original formulation, the gas in question was HI. We now know, of course, that the relevant 
clouds are molecular. 



- 13 - 



We may conveniently recast this fornula in terms of the effective sound speed: 

M, = e-^X . (30) 
Here, the nondimensional quantity 1 is expressed using the traditional, polytropic variables: 

l-2n Ko 

X^^^^in-i/2)^o e-^n+i)^edi. (31) 



47r 







Our prescription gives a finite star formation rate for clouds of arbitrarily low density. 
This is clearly an oversimplification. There is no evidence, for example, that HI clouds 
form stars at all. Even within the molecular domain, it may be that stars form only above 



some threshold density. In their study of the Rosette cloud complex, IWilliams et al.l ( 119951 ) 
found that only clumps which are strongly self-gravitating (as assessed by comparison of 
velocity dispersions, masses, and sizes) have internal stars. While a more complete model 
should account for this threshold effect, we shall not include it explicitly, but simply limit 
our discussion to self-gravitating clumps (/? > /?min) that are capable of forming stars. 



4.2. Comparison with Observations 

In Pa per I, we empirica lly determined the star formation history of the ONC. The 

datab ase of lHillenbrandl (119971 ). together with theoretical pre-main-sequence tracks ( jPalla fc Stahler 
1999), allowed us to assign masses and contraction agesjfl In the detailed analysis, we re- 
stricted our attention to the 244 members with M^, > 0.4 Mq; the Trapezium stars them- 
selves have > 7 Mq, and are thus already on the main sequence. This sample is sta- 
tistically complete, in the sense that the oldest stars do not fall below the observational 
sensitivity limit. 

Based on our results from Paper I, Figure 7 shows M*(t), the mass production per unit 
time, as a function of stellar age. Here, we have binned the data in age intervals of 10® yr. We 
have also extrapolated from our subsample to all stellar masses. We did so by multiplying 
the accumulated mass at each epoch by a factor of 1.3. This factor accounts for th e miss ing 
stars with M* < 0.4 Mq, according to the field-star initial mass function of I Scald (119981 ). 



In order to use equation (30) to describe the ONC, we need the effective sound speed 
a-r as a function of time. Our numerical model gives the nondimensional functional relation 



^Recent ly, 'jefFries' has redetermined the ONC distance as 390 pc, rather than the 470 pc used by 

Hillenbran d (1997). If correct, this distance shift will systematically lower stellar luminosities, and therefore 
increase their ages. 
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a{riT) (recall Fig. 3). Similarly, the quantity I contains ipo and ^o, which we also know 
as functions of r/r. Converting these relations to dimensional form requires that we set the 
cloud mass Mo and background pressure Po- 

We now make the critical assumption that the parent cloud of the ONC was, just prior to 
its dispersal, in the minimum-enthalpy state. Then equations (9) and (10) may be combined 
to yield 

Mo = /i ^ (32) 



and 



Po = h 



G 

GMl 



(33) 



where Ro and qt refer to the final cloud state. The results of our numerical integration give 
/i = 2.0 and /2 = 0.028. We take i?o to equal the radius of the stellar cluster (2.5 pc; see Hil- 
lenbrand 1997), and identify ar with the observed ste/Zar velocity dispersion (2.4 km s~^; 
Jones & Walker 1988). We then find that Mo = 6900 M© and Po = 1-1 x 10"^° dyne cm 
The latter is about 300 times the canonical value in the diffuse interstellar medium, i.e., that 



see 

2 



bounding HI clouds (IWolfire et al.lll995l ). 



It remains only to adjust e, rj, and the exponent n, until the theoretical star formation 
rate M*, as given by equation (30), matches the empirica l one. We use a s tandard imple- 
mentation of the Levenberg-Marquardt fitting algorithm (jPress et al.lll988l . § 15.5). The 
likelihood function that we maximize incorporates a uniform error in the star formation rate 
at each epoch of a = 3 Mq Myr~^. Here we have used a ~ \/N (M), where is the median 
number of stars produced per 10^ yr, and (M) is the average stellar mass at the appropriate 
epoch. This formula assumes that the number of stars in each bin is Poisson distributed 
about the mean predicted by the model, i.e., we neglect the observational contribution to 
the error. 

The dashed curve in Figure 7 shows our theoretical rate as a function of stellar age, 
along with the optimal values for the three parameters. As predicted, e is small (2 x 10~^), 
signifying a low efficiency for stellar production. Specifically, the cloud converts 8 percent of 
its mass (550 Mq) into low-mass stars. The parameter r] is also small (0.3), indicating that 
the cloud contracts over an interval long compared with the free-fall time, tg. The latter is 
1.6 Myr for our initial cloud state. Finally, the best-fit value of n (1.4) lies close to unity. 
Thus, the original star formation law of Schmidt may hold quite generally. The 1-sigma 



^This figure is high er if many of the embedded, near-infrared sources s een behind the ONC arc members 
of the original cluster (|Ali fc Depovlll995t iHillenbrand fc Carpenteiil2000( ). 
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errors in e, rj, and n, are, respectively, 1 x 10"^, 0.08, and 0.09. These fi ffures would have 
been larger had we included observational sources of error. 

We note also that the basic physical characteristics of our model clouds are consistent 
with the clump properties inferred from observations. Our fiducial density po is equivalent 
to a molecular hydrogen number density of 520 cm~^. As can be seen in Figure 5, the 
average interior density is higher by an order of magnitude. The mean visual extinction of a 
cloud, using Mo/n as th e typical colum n dens ity, is 11 mag. These figures are in general 



accord with the findings of IWilliams et al.l (1l995l ) for self-gravitating clumps in the Rosette 



complex. Moreover, Ay, measured inward from the edge, quickly exceeds unity in all our 
models. Thus, the hydrogen is indeed in molecular form throughout the bulk of the interior, 
as is appropriate for a star-forming cloud. 



Discussion 



In this paper we have adopted a conceptually simple model of cloud contraction and 
stellar group formation. The cloud is a self-gravitating sphere, supported against collapse 
by the motion of turbulent eddies. We assigned an effective pressure to this motion; our 
formulation implicitly assumes the eddy speed to be spatially constant, though varying in 
time. Cloud evolution is mediated by the slow leakage of energy, presumed to occur through 
internal, shock dissipation. 

This model, supplemented by a Schmidt-type prescription for the star formation, can 
account not only for the empirically known history of the ONC, but also for more general 
characterist ics of stellar birth. Consid er, for example, the issue of formation efficiency. In 
the spirit of iKrumholz &: McKed (120051 ) . we may define a nondimensional star formation rate 
per free-fall time: 

Cff = TT- ■ (34) 

f M 

Lev i o 

Here, tev is the time over which the cloud produces stars, while AM* is the total mass in 
these objects. Using our ONC model, and setting tg^ = 1 x 10'' yr, we find eg = 0.014. 



Zuckerman fc Evand (Il974l ) long ago pointed out that only about 1 percent of the Galaxy's 
molecular gas can become stars in a cloud free-fall time, to reproduce the observed, global star 
formation rate. The agreement here suggests that giant complexes create stars principally 
through their slowly contracting, internal clumps, as we have modeled. 

We have not described, in any quantitative way, the physics underlying the turbulent 
dissipation. Our best-fit value of rj for the ON C is similar to that found in numerical 
simulations of turbulent clouds (IMac Lowl Il999l ). However, our physical picture is quite 
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different. All simulations to date, which focus on an isolated, interior volume, find the 
turbulent energy dying away. In our model, the mean turbulent speed increases with time 
(see Figure 3). Future, global simlations of self-gravitating clouds supprted by turbulent 
pressure should show this effect. 

Our best-f it value of n agrees, perhaps fortuitously, with that originally proposed by 



Schmidtl (119591 ). For an n- value of unity, the star formation rate per cloud mass scales with 
the gas density. Did this proportionality really hold in the ONC? The righthand panel of 
Figure 5 suggests that it did, at least roughly. Here, the solid curve is the density profile 
of our final, minimum-enthalpy, cloud model. The dot-dashed curve is the current number 
density of ONC stars, as reconstructed from the observed, projected number density (see 
Fig. 3 of Paper I). The similarity of the two curves indeed suggests that stars trace the mass 
distribution of the parent cloud. The same point is evid ent when comparing t he projected 
stellar density with CO contours of the remnant gas (see lStahler fc Fallal l2004l . Fig. 12.27). 



Of course, all stars travel some distance from their birth sites. They do not move 
ballistically, but are subject to the gravitational potential of the parent cloud. Because of 
the star formation law expressed in equation (27), stellar births are indeed concentrated 
toward the cloud center, but there will inevitably be some outward diffusion. In a future 
paper, we hope to track this process through a direct, numerical simulation. 

The minimum-enthalpy cloud that terminates our dynamical sequence is dynamically 
unstable. How do we interpret this instability in a more realistic setting? The essential 
fact is that the self-gravity of the gas becomes so strong that it leads to rapid, internal 
contraction, perhaps even true collapse of the central region. It is tempting to link this 
event with the formation of high-mass stars. While the physics of massive star formation 
is far from clear, the collapse or coalescence of den se, gaseous structures appears to play 
a key role (jStahler et al.l l200d : iMcKee fc Tarul2003l ). In the specific case of the ONC, the 
Trapezium stars are, o f course, centrally located, and appear to be of relatively recent origin 
JPalla fc StahleilboOlh . 



In comparing our model with data from the ONC, we have accepted at face value the 
stellar ages inferred from the placement of each object i n the HR diagram . The age spread 
within clusters remains a contentious issue. For the ONC. lPalla et al.l (120051 ) have found that 
four nominally older stars are de pleted in lithium, as would be expected. Such findings are 
inconsistent with the assertion by iHartmannl (120011 ) that the ostensible age spread primarily 
reflects observational uncertainties. 



Pflamm-Altenburg &: Kroupal (120071 ) accept the higher ages, but hypothesize that all 
such stars were gravitationally captured from somewhat older, neighboring clusters. As the 



-17- 



authors themselves note, the existence of such neighboring systems is unclear. The other 
subassociations within the Orion complex are too young. One possibility is that there were a 
large number of nearby small groups producing low-mass stars and then dispersing. Pending 
more direct evidence for such groups, we continue to believe that the nominal age spreads 
in both the ONC and other systems are real. 

In the present model, we have taken the cloud to be of fixed mass. This assumption may 
be acceptable for the ONC progenitor cloud, at least until the point when the Trapezium 
stars ionized and drove off the gas. In clouds producing low-mass T associations, the latter 
process docs not occur. Yet these clouds arc still dispersed, presumably through stellar 
winds. In our next paper, we will generalize our model of cloud contraction to include the 
effect of continuous mass loss. We will thus achieve a fuller picture of stellar group formation, 
a process of importance not only locally, but on galactic scales. 

We are grateful to Steve Shore for pointing out the work of P. Chavanis on the thermody- 
namics of self-gravitating spheres. This project was supported by NSF grant AST-0639743. 
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Fig. 1. — Evolution of the cloud's internal structure. Shown are the radii of selected, La- 
grangian mass shells as a function of the density contrast (3. From left to right, the five 
shells enclose 0.2, 0.4, 0.6, 0.8, and 1.0 times the total cloud mass. Also indicated are the 
minimum /3-value for self-gravitating clouds and the maximum value for dynamical stability. 
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Fig. 2. — The run of specific enthalpy h along the cloud sequence. This quantity is plotted 
as a function of the effective sound speed a. The dotted portion of the curve pertains to 
clouds which have too low a density contrast to be self-gravitating. 
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Fig. 3. — Evolution of the cloud radius {solid curve) and the efTective sound speed {dashed 
curve). Notice that the time coordinate rjT starts at the first self-gravitating configuration, 
i.e., that for which /3 = Pmm- 
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Fig. 4. — Evolution of the nondimensional central density. As in Figure 3, the time rjr is 
measured from the first self-gravitating cloud. 
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Fig. 5. — Left panel: Evolution of the density as a function of radius. From bottom to top, 
the corresponding values of r]T are 0, 0.60, 0.92, and 0.96. The latter value corresponds to the 
minimum-enthalpy state. Right panel: The density profile of the minimum-enthalpy state 
{solid curve) compared to the reconstructed stellar number density in the ONC {dot- dashed 
curve) . 
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Fig. 6. — Evolution of the velocity profile. In order of deepening minima, the curves corre- 
spond to Tyr-values of 0.60, 0.92, and 0.96. Note that the velocity is normalized to rja. 
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Fig. 7. — Total star formation rate in the ONC as a function of time. The latter is actually 
shown as the stellar age. The solid histogram uses the empirical ages from Paper 1, binned in 
1 Myr intervals. The dashed curve is the theoretical prediction. Also shown are the best-fit 
values for the model's three free parameters: rj, e, and n. 



